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FOREWORD 

Part of this work is taken from a doctoral dissertation submitted by David L. Blanchard, and 
directed by F. K. Chan, to the Department of Space Science, and Applied Physics, Catholic University 
of America, Washington, D.C. The work is an outgrowth of the course “Advanced Celestial 
Mechanics” taught by the second author during the academic year 1968-69. 

Basic notation is as follows: 

(1) V always denotes a column vector and V r denotes the row vector that is the transpose of 
V. 

(2) A r , A -1 , and A + always denote the transpose, inverse, and adjoint, respectively, of the 
matrix A. (Neither B* nor B' denotes the complex conjugate, adjoint, or transpose of the 
matrix B.) 

Diagrams in the “Results” section of this document are from Methods of Theoretical Physics, 
Part I, by Morse and Feshbach, copyright 1953 by McGraw-Hill, Inc., and are used with per- 
mission of the McGraw-Hill Book Co. 
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SEPARATION OF VARIABLES IN THE SPECIAL DIAGONAL HAMILTON-JACOBI 
EQUATION-APPLICATION TO THE DYNAMICAL PROBLEM 
OF A PARTICLE CONSTRAINED ON A MOVING SURFACE 

David L. Blanchard 
Goddard Space Flight Center 

i, and 

F. K. Chan* 

Catholic University of America 


INTRODUCTION 
Historical Background 

The birth of a systematic approach to mechanics essentially occurred with Newton’s statement of 
the fundamental laws of motion in his work Philosophiae Naturalis Principia Mathematica (1687). 
Aristotle (384-322 B.C.) was aware of the principle of the lever but erred in thinking of “virtual 
velocities” rather than “virtual displacements.” Archimedes (287-212 B.C.) was the first to state 
the law of statics correctly; however, in the formulation of dynamics, there was virtually no 
progress for the next 2000 years because all the scholars of mechanics worked with the false 
concepts of the Aristotelian school that associated forces with all motion. It was Galileo (1 564- 
1642) who first became aware of the principle of inertia when he quite accidentally studied motion 
of a body down an inclined plane and out of a horizontal plane. 

Newton (1642-1726) stated his three basic laws as follows: 

( 1 ) 

( 2 ) 


(3) 


Every body continues in its state of rest or uniform motion in a straight line unless it is 
compelled to change that state by forces imposed upon it. 

The rate of change of momentum is proportional to the impressed force and is in the 
direction in which the force acts: 

— = F 
dt 

To every action there is always opposed an equal reaction. 


( 1 ) 


* Present affiliation: Computer Sciences Corporation. 



These three laws have provided the basis for the subsequent development of (Newtonian) mechanics. 
(For an informative account of the history of mechanics, see ref. 1 .) A branch of mechanics developed 
strictly from Newton’s laws of motion is now usually referred to as “vectorial mechanics” and its basic 
viewpoint is to account for all the forces acting on a body and its consequent motion under these 
forces. The basic formulation of mechanics remained unchanged for about 100 years after its 
introduction by Newton. 

It was Lagrange (1736-1813) who gave mechanics a second great forward step when he intro- 
duced the concept of generalized coordinates and a formulation that (under fairly general conditions) 
depends only on a knowledge of the kinetic and potential energy. This branch of mechanics is usually 
called analytical mechanics. Partial credit for its development must also be given to Leibniz (1646- 
1716), Euler (1707-83), and D’Alembert (1717-83). Of course, both analytical mechanics and vectorial 
mechanics result in the same equations of motion. Lagrange’s formulation is independent of physical 
or geometrical considerations once the two scalar quantities, the kinetic and potential energy, are 
known in analytical form. As he so modestly says in the preface of his book , Mecanique Analytique 
(1788), “The reader will find no figures in the work. The methods I set forth do not require either 
constructions or geometrical or mechanical reasonings: but only algebraic operations, subject to a 
regular and uniform rule of procedures.” This was the springboard of analytical mechanics and is sum- 
marized by his differential equations for a system with n degrees of freedom: 

d l ZEX ZE. 

- = Qj i=l,2,...,n (2) 

dt \ZqJ 9 q. 

where 

E k = kinetic energy 
q. = zth generalized coordinate 
Q x = zth component of the generalized force 

About 50 years later (1834-35), the third great jump in the development of mechanics was made by 
William Rowan Hamilton (1805-65). (Some work was also done by Poisson (1781-1840) and Lagrange.) 
By a Legendre transformation he converted (refs. 2 and 3) the set of Lagrange’s second-order differen- 
tial equations (2) into a set of 2 n first-order differential equations. These new equations are called 
Hamilton’s canonical equations and they treat momenta and position on an equivalent basis as the 
independent variables. The well-known Hamilton equations of motion are 

• 9 # 

Qi= r— z = 1 , 2, . . . , n (3a) 

3 Pi 


Pi = 


where 


9 H 

3 q t 


i = 1 , 2, . . . , n 


(3b) 


H = H( q, p, t) = the Hamiltonian 

p i = dL/'dc'jj = generalized momentum associated with the zth generalized coordinate 


2 



L = E k - E p = the Lagrangian 
E p - potential energy 

Next, for the case of a Jime-independent Hamiltonian, by considering the action integral and constant 
energy surfaces H = h in the (q, p) space, Hamilton (refs. 2 and 3) showed that 

Pi = r~ W( q, q') (4a) 

p'i = ~r r^q.q') (4b) 

where (q', p') denotes the initial values of (q, p) at time t = and that the function lV(q, q') satisfies 
the two simultaneous partial differential equations 



- h = 0 


- h = 0 


(5a) 

(5b) 


He called the function W the principal function. For the time-dependent case, by considering the ex- 
tended action integral and the condition that the extended Hamiltonian function H + p t equal zero for 
all t, Hamilton then showed that equations (4) and (5) are to be replaced by 


Pi 


31V 


(6a) 


and 



where IV = lV(q, q', /, t'), which he still referred to as the principal function. Hamilton observed that if 
lV(q, q', t, t') is known, then by simple differentiation and eliminations in equations (6a) and (7a), 
the following equations are obtained 


3 



(9a) 


q i = <7/(q\ p\ t; t ) 

Pi = Pj( q', p\ t'\ t ) (9b) 

(Equations (6b) and (7b) are automatically satisfied by virtue of the condition that H + p ( = 0 for all t 
and equations (8).) The solution of the dynamical problem is then known in terms of the initial con- 
ditions (q', p', t') and time t. This is one of Hamilton’s outstanding discoveries (ref. 4), even though he 
did not give an account of how to determine the principal function. 


Jacobi (1804-51), by considering canonical transformations and generating functions, gave the 
fourth and probably the most gigantic step in the development of mechanics. He showed (ref. 5) in 
1837 that one need only consider one partial differential equation in the generating function S(q, a, t ), 
which is presently known as Hamilton’s principal function 


, bS \ dS 
H [q, — ,t) + — = 0 


0q 


bt 


( 10 ) 


Moreover, he proved that any complete integral of equation ( 10); i.e., a solution with the property that 


d 2 S 
3a. bq. 


? 0 


will provide the solution of Hamilton’s equations (3). (See subsection entitled “Formulation of the 
Hamilton-Jacobi Equation.”) However, Jacobi did not give any systematic rule for finding complete 
integrals of first-order partial differential equations or even, in particular, for equation (10). He did 
seek solutions of equation (10) that could be expressed as the sum of functions each of which is a 
function of only one coordinate; i.e., solutions of the form 


n 

S(q,r, 0 t)= ^ x.{q ,.;«)- T(t;a) (11) 

/= l 

By seeking solutions of such a form, he meant the “separation of variables.” As a tribute to the con- 
tributions of Jacobi, the partial differential equation (10) is now referred to as the “Hamilton-Jacobi” 
equation. 

Much of the subsequent work on the separation of variables of the Hamilton-Jacobi equation in- 
volved the case for which the Hamiltonian is a constant. For this case it is easy to show that there 
exists a solution of the form 


S(q, t;a)= W(q;a y ,a 2 ,. . .,<*„)- (12) 

where the function W is now known as Hamilton’s characteristic function which satisfies the partial 
differential equation 

/ 01V\ 

//(,.-)= a, 03) 

In the period 1891-93, Stackel (refs. 6 and 7) established a theorem giving a necessary and suf- 
ficient condition on the coefficients a ; *(q) and E(q) for the separation of variables for the Hamilton- 
Jacobi equation, if the Hamiltonian is a constant and of the special form 
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( 14 ) 


h = 2 ^(q)p, 2 + nq) 

1=1 

Shortly after this, in 1904, Levi-Civita (in a letter to Stackel; also ref. 8) extended the results to include 
a more general Hamiltonian of the form 

n n 

H = ^ a*{q)p i p j + f^(q) (15) 

/=i /= l 

However, Levi-Civita’s results are not expressed in an explicit form as Stackers are, but are in terms of 
a system of partial differential equations to be satisfied by the coefficients aS(q) and F(q). Conse- 
quently, they are of very little direct practical use. 

For the Hamilton-Jacobi equation, there are infinitely many coordinate systems that have coeffi- 
cients a?(q) that satisfy the conditions obtained by Stackel. Many researchers next used the results of 
Stackel to study the Schrodinger equation. In 1927 Robertson (ref. 9) established a necessary and 
sufficient condition for separation of variables for the Schrodinger equation in the form of a product 
of functions, each of which is a function of only one coordinate. He found that to effect this separa- 
tion, the coefficients in the Schrodinger equation have to satisfy Stackers condition on the coefficients 
in the Hamilton-Jacobi equation and also an additional condition which is usually referred to as Robert- 
son’s condition. In 1934 Eisenhart (ref. 10) proved that there are only 1 1 real coordinate systems that 
are of the Stackel form in Euclidean three-dimensional space; moreover, these 1 1 systems automatically 
satisfy Robertson’s condition. 


Formulation of the Hamiltonian 


To make this work self-contained, the next step will be to set up the Hamiltonian function under 
rather general conditions as a function of n generalized coordinates, the corresponding momenta, and 
time. This explicit form allows one to write out the Hamiltonian function once the kinetic and poten- 
tial energy have been expressed analytically. Furthermore, it is quite easy to relate the Hamiltonian in 
this form to the energy of the system. In the process, the basic principles of mechanics that lead to 
this very general form of the Hamiltonian will be reviewed and nomenclature to be used later will also 
be introduced. A system of N particles will be considered; other systems, such as those consisting of 
many finite parts connected by complicated constraints, can also be included without much difficulty. 

The time-dependent Hamiltonian will be developed with only the following assumptions: 


( 1 ) 

( 2 ) 


The k constraints are holonomic and the N particles have n(= 3N- k ) degrees of freedom. 


The generalized forces Q t are derivable from a scalar work function U through an equation 
of the form 


Q{( q, q, 0 = - 



d /3U\ 
dt \dqj 


(16) 
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where 


n 

U(q, q, /) = E p (q, t) + h,'(q, t)q ( 

«'= 1 


( 17 ) 


Forces described by equations (16) and (17) are said to be monogenic, and it is obvious that conserva- 
tive forces are a special case for which U = E p (q) and b'(q, t) = 0 in equation (17). (In contrast to the 
definition of monogenic forces here, Lanczos (ref. 4) defines them by eq. (16) only. The additional 
restriction given by equation ( 1 7) is included here to eliminate forces dependent on accelerations, 
which do not seem possible in Newtonian mechanics.) 


One could proceed here and use vectorial mechanics and Newton’s laws (eq. (1)); however, when 
the coordinate system is not Cartesian or when constraints are present, Newton’s equations are quite 
cumbersome. It turns out that the approach of analytical mechanics is more suitable. Beginning with 
assumption (1), the equations of transformation to generalized coordinates are given by 

r*=r*(<7 i><7 2 > •••><?„> 0 k = 1,2, . . . ,N (18) 

From equation ( 18), the velocity is obtained by the rules of partial differentiation as follows: 


Kinetic energy is defined by 


d _h 



dt 



n 

S 


!!i + 3r * 

z 

1=1 


dt dt 


1 N 

E k= ; I m k\ * v * 

*=i 


1 ^ dx k dx k 

= ? z m * 

2 dt dt 


(19) 


( 20 ) 


Combining equations (19) and (20) yields the following expression for the kinetic energy in terms 
of the generalized coordinates, generalized velocities, and time: 


n n n 

i=i i=i /= i 

where 

a = a(q, t) 


K 9r * 

— 

dt dt 



( 21 ) 


(22a) 
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a'i = a'.( q, t ) 


3r * 3r * 
> w* — 

* = i 3/ 3<7 f . 


(22b) 


% = V q > 0 


A dr. dr. 


1 v * Wl k 

- — y rn. — • — 

2 k=i M t Mj 


(22c) 


Lagrange’s equations may be derived in several ways, one of which is through D’Alembert’s principle. 
(See ref. 11, 12, or 13.) In one of their more general forms, Lagrange’s equations may be written as 


d bE\ d E k 


dt \d^y d q t 


i= 1,2, . . . ,n 


where the components of the generalized forces are defined by 

JE dr. 

or L*rf 

Fi b( ii 

From equation (16) in assumption (2), it follows that Lagrange’s equations become 


where 


d_ /dL\ _ 

dt \d qj bq. ° 


i = 1, 2, . . . ,n 


L =E k - U 


New variables p., which are called the generalized momenta and are conjugate to the generalized 
coordinates q., are introduced and defined by the equation 


Pi = TT L(q, q, t) 

Mi 

A new function H( q, p, t) is introduced and defined by 


i = 1 , 2, . . . , n 


“ft 

H(q,p,t)= ^ - L(q, q, t) 

_i'=i Jq=q 


L/=i Jq=q(q,P,D 

where q = q(q, p, t) denotes the set of equations obtained from the inversion of the set of equations in 
equation (26). Then, it can be shown (for example, see refs. 1 1 to 13) that the equations of motion are 
given by 

<7.- = ~~H(q, p, t) (28a) 

9 Pi 



(28b) 


d 

Pj = - r— H(q, p, t ) 

Mi 

which are known as Hamilton’s equations. 

From equations (17), (21), and (25), the Lagrangian can be written explicitly as 


L{ q, q, t) = a + 

or, in matrix form, 
where 


Z 2 2 “ipfir E p- 

i=l i'=l /=l 

L = q r Aq + b r q - U 

A 


z b '& 

(29) 

1=1 

(30) 


(31a) 


(31b) 


(31c) 


From equations (22), it is seen that 

(32) 


b = a' - b' 


U = E p -a 

and the superscript T denotes the transpose of a vector or a matrix 

v A A = A 7 


a . = a . or 

// /i u 


From equations (26) and (30) to (32), it follows that 

p = 2Aq + b (33) 

Solving equation (33) for q yields 

> q=-A _1 (p-b) (34) 


where the superscript - 1 denotes the inverse of a matrix. Using equations (30) and (33), equation (27) 
becomes 


//(q,p, /) = (q r Aq+(^. =q(q>p ?) 

(35) 

Then substitution of equations (34) and (32) into equation (35) gives 


H(q, p , 0 = ^ (p - b) r A~ 1 (p - b) + U 

(36a) 

or 


H( q, p, t ) = p T A*p + b* r • p + V 

(36b) 

where by definition 


i , i K/i + 

[a*] = A* = - A -1 = — 

(/ 4 4 IAI 

(37a) 

b*r = - — b r A" 1 
2 

(37b) 
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V= U+ — b r A _1 b 
4 


(37c) 


and the superscript + denotes the adjoint of a matrix. 

Equations (36) are the most general form of the Hamiltonian that are possible under the two 
assumptions stated in this section. One can write the Hamiltonian directly after the equations of trans- 
formation (18) are used to obtain the coefficients a t - given in equations (22) and the scalar work func- 
tion U of equation (17) is specified. In general, this Hamiltonian may contain both quadratic and 
first-order terms in the momenta. It is interesting to note in passing that even if the scalar work func- 
tion is assumed to include second-degree terms in q. t then the Hamiltonian 7/(q, p, t) would still be of 
the form given by equations (36) except that the matrix A would include terms from the scalar work 
function U. 

Finally, a comparison is made of the Hamiltonian H and the total energy h defined to be the sum 
of the kinetic energy and the potential energy. To discuss meaningfully the usual concept of potential 
energy, it is noted from equations (16) and (17) that b' must equal zero and the potential energy can 
then be taken to be E . Consequently, from the expression of the kinetic energy E k given by equation 
( 21 ) 

n n n 

h= 2 S W/ + S a 'A +a + E P (38) 

i= 1 /= 1 /= 1 

Substitution of equations (29) and (33) into (27) results in 

I«q,P.O] p . p „ i4> „= £ £vW/-« + £, i> <39) 

1=1 /= 1 

which remains valid regardless of any requirement on b'. Hence, from these equations, it follows that 
the total energy is identical with the Hamiltonian if and only if 


o' i q i + 2a = 0 (40) 

i= 1 


where a and a' are defined by equations (22a) and (22b). It is obvious that a sufficient condition is 
given by 


a = 


a 


1 

i 


1 f 

3r * 

9r * 


m. — • 

— = 0 

2 La 

k= 1 

* at 

dt 


JL 3r. 0r 

> m, ~ " = 0 

IR 3 9 1 9<7,. 


(41a) 


(41b) 


These conditions are satisfied when the equations of transformation (18) are time independent; i.e., 
for coordinate systems at rest. 


9 



Formulation of the Hamilton-Jacobi Equation 

Having obtained Hamilton’s equations 


<7/ 


3 

—H( q, p, 0 
9 Pi 


3 

Pi = ~ #(q, P> 0 

the next step is to introduce new coordinates and momenta (Q, P) defined by 

Q = Q(q, p, t) 


(28a) 

(28b) 


(42a) 


P = P(q, p, t ) 


(42b) 


These new variables must be such that there exists a function K( Q, P, t) for which the new equations of 
motion are given by 

4=^-m p ,0 (43a) 

^ ' 

p i = - p ’ '> (43b) 

oQ t 

where AfQ, P, t) is known as the new Hamiltonian. Transformations (42) for which equations (43) 
hold are known as canonical transformations. (In older literature the term “contact transformation” de- 
notes the case for which time does not appear explicitly in equations (42) and the term “canonical 
transformation” denotes the case for which it does. Such unnecessary differentiation will not be made 
here.) Equations (28) constitute the set of Euler equations corresponding to the condition 




pfii - H(q, p, t) 


dt = 0 


whereas equations (43) correspond to 


(44) 


5 




- K(Q, P, 0 


dt = 0 


(45) 


Hence, given equation (44) or (28), it is sufficient (but not necessary) that there exists an arbitrary 
function F such that 


2 Pfli - P. = ^ P £i - P. 0 + — (46) 

i^l i= 1 


in order for equation (45) or (43) to hold. That is, transformations (42) for which equation (46) holds 
are canonical, but not all canonical transformations are given by equation (46). The fact that equation 
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(46) is not a necessary condition is easily seen by noting that given equation (44), then transforma- 
tions (42) for which 

- K(Q, P, t) (47) 

where c is an arbitrary constant, are also canonical. The problem of finding a necessary and sufficient 
condition for canonical transformations as defined by equations (43) is not important to the solution 
of dynamical problems. All that is needed is a sufficient condition such as that given by equation (46). 
(There is some confusion in the literature about equation (46) being a sufficient or necessary condition 
for canonical transformations. For instance, Goldstein (ref. 1 1) erroneously states that it is necessary, 
whereas Nordheim and Fues (ref. 14) and Born (ref. 15) state that it is necessary and sufficient. 
Gelfand and Fomin (ref. 1 6) state it correctly and so does Lanczos (ref. 4) except that at the end of 
his discussion he calls it the most general condition for a canonical transformation, and thus causes 
undue confusion. Another important point of confusion is that some authors, e.g., Pars (ref. 12) 
and Whittaker (ref. 13), define canonical transformations for which equation (46) is valid and 
equations (43) are not. In this case, there is no question to begin with as to whether equation 
(46) is a necessary or a sufficient condition for a canonical transformation.) The function F 
appearing in equation (46) is called the generating function. To effect a transformation between the 
new variables (Q, P, t ) and the old variables (q, p, t), F must be a function of one of the following 
forms (see ref. 11 or 1 5): 


F = F l (q, Q, t) 

(48a) 

F = F 2 (q, P, t) 

(48b) 

F= F 3 ( p, Q, t) 

(48c) 


Z *W-ff(q,p,/) = C y P t Qf 

i=l [m 


F = F 4 (p,P,0 (48d) 

In any of these cases, the new Hamiltonian K{ Q, P, t ) and the old H( q, p, t ) are related by 

9F 

K = H+ — (49) 

dt 

The next step is determination of canonical transformations such that the new Hamiltonian K is a 
constant, which without any loss of generality can be taken to be zero. Therefore, it follows that the 
generating function F must satisfy the equation 


9 F 


H(<\, p, t) + — = 0 
01 


(50) 


For all mechanical systems H (as given by equation (36)) is known in explicit form for the momenta 
p but not for the coordinates q. Consequently, either the transformation of equation (48a) or (48b) 
is the best choice of F. For equation (48a) 
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(51a) 


P. = -■ 

l 


(51b) 


and for equation (48b) 


Qi= ~ 
a p : 


(52a) 


(52b) 


In either case, by combining equations (50) and (51) or (50) and (52), it follows that Fmust satisfy the 
partial differential equation 

/ dF \ dF 
H[q, ~-,t) + — =0 


If F is replaced by S so as to conform to common usage, the equation becomes 

„ / as \ as 

H (q. + TT =0 


which is known as the Hamilton-Jacobi equation. Some authors prefer to use S = F 2 (q, P, t ) and 
others S = F l (q, Q, t ), but, as previously demonstrated, there is no real difference between these 
choices. In view of equation (36b), the most general Hamilton-Jacobi equation (53) may also be 
written explicitly as 

v v as as ^ as as 

Z 2 a ^ (q ’ 0 a “^ + 2 6;(q ’ 0 + nq ’ 0 + ^7 = 0 


i= 1 /= 1 


Jacobi (ref. 5) showed that the solution of the dynamical problem defined by Hamilton’s equa- 
tions (28) can also be obtained by finding any one complete integral, denoted by S = S(q, a, t), of 
equation (53). 

Using equation (48a) and letting new coordinates Q be the arbitrary constants a and the new 
momenta P be the constants /3, one obtains from equations (51) 

3 

P t = — S(q, a; t) (55 

d Qi 


0/ = - — S(q,a; t) 
3a ( . 


(55b) 
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Equation (47) can then be solved to yield 


ft = ft(«, ft t ) 


(56a) 


Pi = P,(«, ft t) (56b) 

The constants a and (i are obtained from the initial values (q\ p') at time t', but they are not necessar- 
ily identical with the initial values. This is the main difference between Hamilton’s partial differential 
equations (8) and Jacobi’s partial differential equation (10). It is also noted that to invert equation 
(55b) to obtain equation (56a), the following property of complete integrals is a necessary and suffi- 
cient condition: 

a 2 s 


3ft bq j 


^0 


(57) 


Similarly, using equation (48b) and letting P = a and Q = ft from equations (52) one obtains 


Pi = ~S( q, a, t) 
9ft- 


(58a) 


ft - Z S(q, a, t) (58b) 

' da. 

i 

which can be solved to yield a set of equations of the form of equations (56). Hence, the next step is 
to find complete integrals of the Hamilton-Jacobi equation. It follows from this discussion that once 
this is done the evolution of the whole dynamical system will be known. 


Outline of Present Work 

In the next section, analysis is made of the time-dependent Hamilton-Jacobi equation: 

v /a s\ 2 as 

X +K(q ’ 0+ i7 =0 (59) 

for which necessary and sufficient conditions on the coefficients ft*(q, 0 and F(q, t) are sought so that 
there exists a complete integral of the form 

n 

S(q, t; a) = £ Xfa ; a) - T(f, a) (11) 

Because equation (59) does not contain cross terms (dS/dq^ibS/dq.) or first-degree terms (as/3ft), for 
convenience it will be referred to as the special diagonal Hamilton-Jacobi equation in contrast to 

/a s\ 2 sr- as as 

< (q - '> y + 2 6 * (q - » 37, + () + 7 = 0 <60) 
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and 


V' dS dS dS 

1 1 w; I + 


ds 

+ — =° 

dt 


(54) 


which are referred to as the general diagonal Hamilton-Jacobi equation and the general Hamilton-Jacobi 
equation, respectively. 

In the section on solutions, the results of this analysis with n = 2 will be applied to the solution of 
a comparatively wide class of dynamical problems involving a particle moving in a Euclidean three- 
dimensional space under the action of external forces but constrained on a moving surface. 


SEPARATION OF THE SPECIAL DIAGONAL HAMILTON-JACOBI EQUATION 

In this section, necessary and sufficient conditions will be established for the separation of the 
special diagonal Hamilton-Jacobi equation in the form of a sum of functions, each of which is a func- 
tion of only one variable together with the arbitrary constants a. The results are a generalization of 
those obtained by Stackel (refs. 6 and 7) because the Hamiltonian is now time dependent. The 
case of two degrees of freedom will be considered specifically before the general case of n degrees 
of freedom. 


Theorem on Separation of Variables for Two Degrees of Freedom 

Consider the equation 

(dS \ 2 /dS \ 2 dS 

a * i(x,r,/) \teJ +a *2 (x ’ y ’ t) \fy) + V( *>y> t)+ j t =0 

where, for convenience, the following notation has been introduced: 

Q \ -x 

q 2 = y 

It is desired to investigate the existence of a complete integral of the form 

S(x,y, X(x-,a l ,a 2 )+ Yiy;^, a 2 ) - T(t;a x ,a 2 ) 

i.e.; a solution with the property that 


d 2 S 

d 2 S 

da j 9x 

dttj dy 

d 2 S 

d 2 S 

da 2 dx 

0a 2 dy 


( 61 ) 

(62a) 

(62b) 

(63) 

(64) 
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where X(x; ctj , a 2 ), Y(y ; , a 2 ), and T(t; ojj , a 2 ) are arbitrary functions of x, y, and t, respectively, 

and ofj and a 2 are arbitrary constants. Substituting equation (63) into (61) yields 



(65) 


For convenience, the symbols 0, <j>, and \p are introduced and defined by 


( dX \ 2 

6(x;a a )= — 

1 * \dx ) 

(p(y-a 1 ,a 2 ) = 



(66a) 

(66b) 


^(/;a 1 ,a 2 ) = 


dT 

dt 


(66c) 


Then equation (65) becomes ' 

a*j0 +a| 2 0 + F= \1/ (67) 

A partial differentiation of this equation with respect to and « 2 yields 

a *l Q \ +a 22^l = (68a) 

a ll 0 2 +fl 22 <S& 2 = ^2 (68b) 

where the subscripts 1 and 2 on 6, 0, and i p only refer to partial derivatives with respect to a, and a 2 , 
respectively. Next, from equations (63) and (66), 


a 2 s 

a 2 s 


a 


3 /dY\ 

a aj ax 

aoij a y 


a«j 

\dx) 

aaj \ dy J 

d 2 S 

a 2 s 


a 


a / dY\ 

da 2 dx 

da 2 3 y 



\dx) 

da 2 \dy ) 



(69) 
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Hence, in view of property (64), it follows that 

0 i *1 

P 0 (70) 

e 2 <p 2 

Consequently, equations (68) may be solved for and a| 2 so as t° obtain 

'Pi (t; a t , a 2 ) 0 2 O; , a 2 ) - 0 2 (t;a t , a 2 )<t> i (y ; a, , a 2 ) 

a* l (x,y,t)= (71a) 

0 2 (x;a 1 ,O! 2 )0 1 (y;a 1 ,Of 2 ) 

- i//j (f; , a 2 )d 2 (x; ^ , a 2 ) + ^ 2 (/; , « 2 )0j (x; , a 2 ) 

tf*2(*>y> 0 = (71b) 

0 1 (jc;ot 1 ,a 2 ) 0 2 O;a 1 ,a 2 )- 0 2 (x;aj , a^O; , a 2 ) 

Because these relations must hold for all values of ccj and a 2 , it follows that a* 2 and a* 2 must have the 
forms 

T[(t)Y 2 (y)~ T 2 (t)Y , (y) 

0 = (72a) 

X { (x)Y 2 (y)- X^Y^) 


~T i (t)X 2 (x) + T 2 (t)X { (x) 

af 2 (x,y, t) = (72b) 

x 1 (x)r 2 (y)-^r 2 (x)y 1 o) 

where X l (x), X 2 (x), TjCy), Y 2 (y), T^t), and T 2 (t) are arbitrary functions of their arguments. 

Next, note that equation (65) may also be written as 

V(x,y, t)= 0(t; o-j , a 2 ) - a* y {x,y, t)0(x; , a 2 ) - a* 2 (x,y, t)<p(y;a t ,a 2 ) (73) 

Again, because this relation must hold for all values of and a 2 , it follows that V must have the form 

V(x,y, t) = T 0 (t)~ a* x (x,y, t)X 0 (x)~ af 2 (x,y, t)Y Q (y) (74) 

where X Q (x), K 0 (y), and T Q (t) are also arbitrary functions of their arguments. 

Finally, it is seen that equations (72) and (74) constitute necessary conditions on a * 1 , a| 2 , and V 
for equation (61) to have a complete integral of the form of equation (63). 

Next, it is shown that these conditions are also sufficient. That is, if a * x , a 22 , and V are of the 
forms of equations (72) and (74), then equation (61) has a solution of the form of equation (63). To 
do this, first consider the equation 


where 


(a 

<X\ 2 (dY \ 2 

dT 


U » Kc 

— + — + V- 

IxJ 22 \dy) 


(65) 


X = X(x;a l ,a 2 ) 


(75a) 
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(75b) 


Y= Y(y;a x ,a 2 ) 


T= T(t;a x , ot 2 ) 

Substituting a* { , a* 2 , and V from equations (72) and (74) into (65) yields 


T, ~ Tj Y, 


(-) - X 

(J ° 


Next, multiply by 


T l X 2 ~ Vl 
X 1 Y 2 ~ X 2 Y \ 


X 1 Y 2~ X 2*X 


fd I ) 2 

\dyj 


dT 


dt 


The result is 
1 

?X X 2 ~ Vi 


(d_xy 

\dx) 


X n 


{T l Y 2 -T 2 Y x )(T x X 2 - T 2 X x ) 

1 


T Y - T Y 

1 1 1 2 I 2 I l 


a- 


It is noted that 


x i Y 2~ X 2 Y i 


(X x Y 2 - x 2 Y x )(dT/dt - T 0 )! 
{T l Y 2 -T 2 Y l ){T x X 2 -T 2 X l ) 

X JT x Y x /T x 


(T, Y 2 - T 2 Y x ){T x X 2 - T 2 X x ) t x x 2 - t 2 x x t x y 2 - t 2 y x 

so that equation (77) becomes 


t x x 2 - t 2 x x |i, 


(l 

r/dx\> ■ 

■ t dT \\ 

1*1 

1 

Tj 

i 



(75c) 


(76) 


(77) 


(78) 


T \ Y 2 ~ T 2 Y x 


i 

K> 

1 

1 fdT \ 

— 

— ) - ^ 

— -T 0 ) 

W 

L \dy/ \ 

T x \dt / 


(79) 


The left-hand side of this equation is a function of x and t only whereas the right-hand side is a 
function of y and t only. This can only be true if each side is equal to the same function of t because 
x, y, and t are all independent variables. 


Consequently, from equation (79), 


T x X 2 - T 2 X x 


’ 1 

~/ dX \ 2 

— 

(- ~ x o 


\dx/ 


T x \dt 


=m 


(80a) 
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T\ Y 2 - T 2 Y X 


q’-r, 

\dyj 


T x W 


=m 


(80b) 


On subtracting these two equations, one obtains 


1 

7 dX\ 2 

1 

’/ dY \ 2 

— 


- — 

( ) — ^0 

*1 

\dx / 

Y i 

\dy) 


= f(t)T l 


X 2 Y x ~ Vx 


XxYx 


or 


r> - Vi 


<-) - x 

U 0 


x , 


X 2 Y X - Vl L 


Am 2 

\dy) 


= f(t)T 1 


(81) 


(82) 


In equation (82), the left-hand side is a function of x and y only, whereas the right-hand side is a func- 
tion of t only. This can only be true if each side is equal to the same constant; therefore, 


and 


/(/)7\ =<*2 


1 

% 

to 

1 

1 

1 

— 

— - *0 

- ' — 

— - ^0 


l\dx/ 

y i 

L \dyf \ 


<x 2 (X 2 Y x - y 2*i> 


(83) 

(84) 


ZiTi 


From equation (84), 

1 


(dX\> 

W " X ° 


X „ 


1 


a. 


fdY > 

Wy/ 


- 


«2 


(85) 


- 1 

Again, the left-hand side is a function of x only and the right-hand side is a function of y only; there- 
fore, they must both equal the same constant. Consequently, 


— ) - X n - a.X 7 = a. X. 
dx) 0 2 2 1 1 


(86a) 


(dY\ 2 

y -iw.-.'. 

Finally, substituting equations (83) and (86a) into equation (80) yields 


X, 


which reduces to 


T X X 2 ~ T 2 X X L 


■ x 2 1 (dT \ 

a 2 + a i " ( T 0j 

X x T l \dt ) 


a. 


dT 

T n ~ ol.T. = a. T 7 

dt 0 112 2 


(86b) 


(87) 


(88) 


Hence, it is seen that if , a| 2 , and V are given by equations (72) and (74), then equation (65) has a 
solution for which the functions X, Y, and T are given by 
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(89a) 



dT 

dt 


= X Q (*) + a 1 X x ix) + a 2 X 2 (x) 

= ^ O 0) + a 1 y i (y) + a 2 r 2^) 

= T 0 (t ) + a 1 T l (t) + a 2 T 2 {t) 


(89b) 

(89c) 


Consequently, it follows that if af 1 , a| 2 , and V are given by equations (72) and (74), then equation 
(61 ) has a solution of the form of equation (63). 

The preceding results are summarized in the following theorem. 

Theorem 1: The equation 


k , (dS \ 2 (bs \ 2 bS 

a*i(x,y, t) [^—J +a* 2 (x,y,t)y—J + V(x,y, t) + — = 0 

has a complete integral of the form 

S(x,y, t;o; 1 ,a 2 ) = X(x;o; 1 ,a 2 )+ T(y;aj,a 2 )- a 2 ) 

rr 

if and only if a*, , a* 2 , and V have the forms 1 

T l {f)Y 2 (y) - ^(OTjO) 

a* n (x,y,t)= ; 

X x {x)Y 2 (y)~ X 2 {x)Y x (y) 

\ 

-T x {t)X 2 {x)^T 2 {t)X x {x) 

a* 2 0 = 

X l (x)Y 2 (y)~ X 2 (x)Y x (y) 


(61) 


(63) 


(72a) 


(72b) 


V(x,y, t)=T 0 (t)~ a* x (x,y, t)X Q (x)- a* 2 (x,y, t)Y 0 (y) (74) 

where X.(x), Y.(y), and T { (t) for i = 0, 1,2 are arbitrary functions of their arguments. The functions 
X(x; aj , a 2 ), T(y; cXj , a 2 ), and T(t; a y ,a 2 ) are obtained from the equations 



= X 0 (x) + tXj X x ix) + a 2 X 2 ix) 
= r 0 (y) + a 1 T 1 O) + a 2 y 2 (y) 


(89a) 

(89b) 


dT 


T 0 ^ + °1 ^1 W + a 2 ^2 W 


(89c) 
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so that the complete integral is given by 


S(x, y, t; aj, a 2 ) = J (X Q + a l X l + a 2 X 2 ) 112 dx 

+ J (7 0 + Oj Kj + a 2 Y 2 ) 112 dy- J (T q + a 2 T x + <* 2 T 2 ) dt 

Theorem on Separation of Variables for n Degrees of Freedom 


In this section the equation 


V /bS\ 2 

1 '» fej 

i= l x 


+ V(q, t) + — =0 


is considered and an investigation is made of the existence of a complete integral of the form 


n 

S(q, t;a)= £ X^ot)- T(t;< 


i.e., a solution with the property that 


-v \ -r v W ' > 

I da / 9< ?/| 

where X i (q i ; a) and T(t; a) are arbitrary functions of q. and t, respectively, and also of 

a= . . . ,aj (91) 

(As a rule, the vector A will always denote a column vector even though it may be written horizontally 
as above for the sake of convenience and conservation of space.) 

Substituting equation (1 1) into (59) yields 

- /dXV dT 

\ a * l +v = 0 (92) 

f=l XdqJ dt 

For convenience, the symbols 0., < p, 9, and 7 are introduced and defined: 


i= \ , 2, ... ,n 


9 i (q i -,a)= 


\p{t\ a) = 


8 = [Sj^jja), 9 2 (q 2 ;a), . . . ,d n (.q n ;oi)\ 
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7 = [afi (q, 0, a* 2 (q, a* n (q, /)] 


Then, equation (92) becomes 


n 

^ 7 fii + V= 4> 


or 


i=i 

7 r • d + V = \}/ 

where the superscript T denotes the transpose of a vector (or a matrix). 

Differentiation of this equation with respect to a f for i = 1,2 , ,n yields 

dO f 3i b 


V / 

> y,- — = 

/=l 3 a. 


/= 1,2, . . . ,n 


or 


where by definition 


0 • 7 = 

0 = %%’> a) 

* = [ (t; a), * 2 (f ; a), . . . , (/; a)] 

"3i|/(t;a) dtp(t;a) di(/(t;a )' 


3a 


3a, 


3a 


From equation ( 1 1 ) it follows that 

a /dx\ 


d 2 S 


da. a q. 


3a, \dq , 


9«, W/ 


a At*, 


3q >, v?l> 


3 “„ Wj 


2 n n — 

i=i 


a 

/dxy 

a«j 

\ dq J 

a 

(dX\ 2 

3a 

n 

\ dq J 


a 

3a, 


fd^ 2 1 


a 

( d K 

da 

n 

\ dq nJ 


3a, 


» dX. 

2 n n ~7~ 

» =l dq ( 


lei 


2 " n 

i=l «<7/ 


(93d) 

(94a) 

(94b) 

(95) 

(96) 

(97a) 

(97b) 


(98) 
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(99) 


Hence, from equations (57) and (98), it is seen that 

101 ? 0 


Solving equation (96) for y yields 


7 = 0“ l • 'if 


Denoting the minor of 0 f . . by . allows equation ( 1 00) to be written as 


7 y (q, t) = • 


n. 

^ (-iy + '* f (t;a)<^. 

1=1 


n 

2 (-1 ,(«/;»)*</ 

i=i 


( 100 ) 


2 (-i r'*, ,(««) ^(±> n 

«=1 P k. Z=1 


P Jfc, /=1 
k*i,l+f 


n 

(-iy' +/ 0 l7 (? /; a) ^ (±) n ®*/(9/;«) 

/> Ar, /= 1 

k¥>i,l¥=j 


i=l 


where for convenience the following notation has been introduced: 
n 


(101) 


£ (1) fl ao2) 


P A,/=l 

k ¥= i, l^j 


in which the summation convention is used, 5* 1 j 1 is the generalized Kronecker delta, and [w^] 


is the matrix of elements appearing in the minor i.e., 


["h! = 






%* di («i;«) 


®„i 


©, w _ 




®<*+lX/-l)ty'-l ; “1 




®l (/> 1)*7/+ 1 ’“1 




®ff-l)(/+l)^/+l ’“1 ' ’ ’ ®<(- linin’®) 


®(/+ 1)0+ 1)^/'+ 1 ’ “1 ' ' ' ®(i+lVi*7/i ’“1 


®n(/+ ljtfyt- 1 > “1 


(103) 
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Because the relations hi equation (101) must hold for all values of a, it follows that 7 (q, t) must have 
the form 


or 


n n 

\ j (±> n x t Mi 


P k, 1=1 

k±i, l*j 


^C-lP'W,) j'w n X t ,(q,) 


P k, 1= 1 

k*i, l*i 

j= 1,2 n 


= X 1 • T 


|X| = ^0 


(104) 


(105) 




(106) 

where 

X = [X.fa)] i,j= l,2,...,n 

(107a) 

\ 


T=[T l (t),T 2 (t),...,T n (t)] 

(107b) 


Consequently, the right-hand side of equation (101) must have numerator ./V and denominator l@l of 
the form 

n n 

N=Mot) £ (-ly+iT^t) 2 (±) FI X kM /) (108b) 

1=1 P k, 1=1 

k*i, /#/ 

101= X(flO J (-ly^.^.) J (±) f[ **,(?,) = X(a) • IXl (108a) 

i=l P k, 1=1 

k*i, l ¥=j 

where X(a) is some arbitrary function that does not vanish identically. 

Equation (94a) may also be written as 

n 

V(q ,/)=*(/;«)- J 7 ,(q, /)0,(^;«) (109) 

i=l 

Again, because this relation must hold for all values of a, it follows that F(q, t) must be of the form 

n 

F(q, t)=T 0 (t)~ ^ 7 ,(q> t)X 0i {q.) 0*0) 

i=i 
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or 


where 


v= T o ~ xT o-y 

*0 - ^01 ^1 X 02^2^ > • • ■ ’ Xonton ^ 


dll) 

(112a) 
(1 12b) 


T Q = T Q it) 

Finally, it is seen that equations (104) and (109) constitute necessary conditions on y and V for 
equation (59) to have a complete integral of the form of equation (1 1). 

Next, to show that these conditions are also sufficient, it is noted that if y is given to have the 
form of equation (105) where X^iq.) and T(f) are arbitrary functions, then &..(q.;oi) and ^.(f; a) may 
be chosen such that 




i, j = 1 , 2, . . . , n 


-w J 


— < p(t; a) = H'.(t;a) 
da. 1 




= T t {t) 

noting that l0l = X(a) • IXl is satisfied with A(a) = 1 . From this, it follows that 


n 

i= 1 


/ = 1 , 2, . . . , n 


(113a) 


(113b) 


(114a) 


n 

Ht; a) = T*(t ) + ^ T f (t)a. 


(114b) 


where X%.(q.) and T'^(t) are arbitrary functions. 

Moreover, if V is given to have the form of equation (110) where X^q^ and T Q (t) are arbitrary 
functions, then from equations (94a) and (101 ) it is found that 

n r n n n n 

y i r v + y w,*, + r 0 (»- y r,jr 0 j < 9 ( )=r*«)+ y tytx, (H 5 ) 

z=l L /-l J i'=l 7=1 

Because 7 (q, t) is arbitrary (to the extent that X j .(q / .) and T.(t) are arbitrary), it follows that 

x l &(*,) = *«<*,) ( 116a > 

T*(t)=T 0 (t) (116b) 
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In passing, for the sake of completeness of discussion, it is noted that because a are arbitrary constants, 
the remaining terms in equation (115) lead to the equation 


1=1 

or 

X • 7 = T 
or 

7 = X 1 * T 

which is automatically satisfied because 7 is given by equation (104) or (105). 
Next, from equations (1 14) and (1 16) it follows that 

n 

d i (q j ; a) = X^) + ^ X^q^. i = 1 , 2, . . . , n 

/= 1 


( 117 ) 

(118) 
(119) 


(120a) 


n 

Mt;a) = T 0 (t) + I 

i= 1 


(120b) 


which are expressions for d i (q.; a) and \jj ( t ; a) that satisfy equation (94). 

Hence, in view of equations (94) and (120), it is seen that equation (59) has a solution of the 
form of equation (1 1) for which X.(q.; a) and T(t; a) are given by 



n 


i = 1 , 2 , . . . , n 


(121a) 


/= 1 


d T 
dt 


r„( 0 + ^ 7)(ooj 

7=1 


(121b) 


where IXl^ 0. The fact that this solution is a complete integral is seen from equation (98) from which 
is obtained 


d 2 S 


101 


3 «- 



dX t 

dq t 


1X1 

=£ 0 

n dX. 

7=1 dq { 


( 122 ) 
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The preceding results are summarized in the following theorem : 
Theorem 2: The equation 

ST (dS\ 2 dS 

I a " (q ’ ,) \^J + F(q ' ,) + v , =0 

has a complete integral of the form 

n 

S(q,t\a)= X j (q i ,a)- T(t,oc) 

i=i 

if and only if y = [a* , a| 2 , . . . , a* n ] and V have the forms 

7 = X-! T 

IXl^O 


( 59 ) 


(ID 


(105) 

(106) 


V= T Q (t)~ Xj • X -1 • T (123) 

where T Q (t) is an arbitrary function and X Q , T, and X are arbitrary vectors or matrix given by 

X 0 = [X Ql 0 ? , ), X 02 (q 2 ), ■ ■ ■ , X 0 n tin O 1 2a > 

T= [Tj (/), T 2 (t T n (/)] (108b) 

X=[X. j (q.)] i,j=\,2,...,n (108a) 

and the functions X.(q.; a) and T(t\ a) are obtained from the equations 

/dX\ 2 

/ = 1, 2, . . . , n (121a) 

(121b) 


— =-w + 2 

\ dc li/ 7=1 

n 

= T 0 (t)+ ^ T.(t)a. 


dT 

dt 


so that the complete integral is given by 



nr n n 

^ ^ 

1/2 r r \r 

^(q, t;a)= ) / 

x «to t )+ 2 W 0 / 

dq t - T Q (t) + yT.(t)a. 

i=lJ 

7 = 1 

J L m J 


(124) 


Remarks 

In the transition from equation (101 ) to (104) it is not required that the arbitrary functions 
0 i/ .(<? / .; a) and T.(t; a) of and t be equal correspondingly to X j .(q.) and T^t). The reason for this is 
that it is possible for additive terms or factors involving a to be eliminated respectively when a sum or 
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a product (or quotient) is considered on the right-hand side of equation (101). The error made by 
Stackel when considering the time-independent Hamilton-Jacobi equation was in stating that 0..(o.; a) 
must be equal to X^{q.). 

It is seen that these results reduce, as one would expect, to those obtained by Stackel for the time- 

independent Hamilton-Jacobi equation merely by letting T x (t)= 1, T ( {t) = 0 for / = 0, 2, 3 n, and 

= h. 


It is possible for the arbitrary constants to appear in a more complicated form in the complete 
integral instead of just the simple form. This is accomplished by taking new arbitrary constants a* 
where 

a* = [a*(a), a*(a), . . . , a*(a)] (125) 

such that 


9a* 


¥ 0 


(126) 


In this case, it is obvious that, in general, © i/ .(<? / .; a) and ^.(f; a) are not equal to X^iq^) and T.(t) and 
that A(a*) in equation (107) is explicitly dependent on a*. 


However, these additional insights do not alter the necessary and sufficient conditions for equa- 
tion (59) to have a complete integral of the form of equation (1 1). Neither do they introduce any 
essential generalization to the complete integral given in theorem 2 where it is stated in its simplest 
form. 


SOLUTIONS OF THE DYNAMICAL PROBLEM OF A PARTICLE CONSTRAINED ON A 
MOVING SURFACE 

In the previous section, conditions were derived for separating the solution of the ^-dimensional 
special diagonal time-dependent Hamilton-Jacobi equation in the form of a sum of functions. It is 
obvious that these results can be applied to the solution of problems in dynamics as long as the con- 
ditions of separation are satisfied. In this section, the comparatively wide class of problems involving 
a particle moving under the action of external forces in a Euclidean three-dimensional space, but con- 
strained on a surface moving as a function of time, will be considered, For these problems, there are 
two generalized coordinates and the coefficients in the Hamilton-Jacobi equation can involve time 
explicitly. Hence, the results stated in theorem 1 may be used. 

The Motion of a Particle Under the Action of Forces in Euclidean Three-Dimensional Space and 
Constrained on a Moving Surface 

Consider a particle moving under the action of forces in a Euclidean three-dimensional space E 
such that its coordinates always satisfy a constraint given by 

f(x,y,z,t) = 0 (127) 

This means that at the instant t = t Q , it lies on a surface 2 whose equation is /(x, y, z, t Q ) = 0; for 
various instants of time 2 will assume different shapes and positions in E. For simplicity, it is assumed 
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that the continuous motion of 2 in E generates a family of smooth surfaces and that no member inter- 
sects with any other member of the family. (See fig. 1 .) One way of relaxing the condition on the 
smoothness of 2 is by considering 2 to be a cone whose angle changes so as to generate a 
family of cones with a common axis and a common vertex. (See subsection entitled “Spherical 
Coordinates” in the “Results” section.) Two ways of relaxing the condition on nonintersecting mem- 
bers of the family are by making 2 be a plane rotating about an axis (see subsection entitled “Cylindri- 
cal Coordinates” in the “Results” section) or by requiring the existence of finite intervals (tj , t 2 ), 

( t 2 , t 3 ), etc., such that no member intersects any other member of the family in the same time interval. 
(See the section entitled “Results” with/(t) a piecewise-monotonic continuous function.) 

Let | and t? be constants defining any two families of intersecting curves always lying on 2. As 2 
moves in E, these two sets of curves will also, in general, change their shape. Thus, the motion of the 
particle in E can be described by the equations 


x = x(X,r},t) (128a) 

y = y(£,r),t) (128b) 

z = zQ, n, t) (128c) 

which can also be interpreted as the equations of transformation from (x,y, z) to generalized coordi- 
nates (£, tj) as a result of eliminating the constraint of equation (127). 

From the introduction it is known that the most general Hamilton-Jacobi equation is given by 



V dS dS 
> a* — — + 

L " a q t 


n 

1 

i=i 


ds as 

&,*— + v+ — =o 

dt 


MOTION OF SURFACE 1 
IN SPACE E 



(54) 


Figure 1.— Geometry of the particle. 
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where a?-, b * , and V are defined by equations (37). Hence, for the problem of the particle constrained 
on a surface moving in the space E, 


Hi 


a S 

3? 


+ a 22 


bS 3 S bS bS bS 

+ 2a* + b* — +b* — + V + — = 0 

12 b V 1 2 brj bt 


(129) 


For simphcity, it is assumed that b' = 0 in equations (16) and (17) so that the external force act- 
ing on the particle is derivable from a potential E . Then, from equations (18), (26), (31), and (128) 
it is seen that 


Ti 


= ^ x t +y ' +z P 


(130a) 


a 22^ j( x v +y i +z i } (130b) 

a l 2 = a 2l = \^ x i x n + y i y n + z t z v ) (130c) 

b x =a\ - b\ =a[ -x ( x ( + y t y^ (130d) 

h = “'i ~ b 'l ‘ a \ ‘ x , x , + + -V% d30e) 


U-E p -a = E p - jOrJ+^+r?) 
From equations (37) and (129) it is found that 


z* = 
'll 


‘11 


4 lAl 


(l/2)(xj+y 2 +zj) 
4 lAl 


I — 

l 22 


‘22 


4 1 Al 


(l/2)(x| +y 2 +z 2 ) 
4 1 Al 


a* = a* — 
a 12 “21 


12 


4 1 Al 


-( 1 / 2 ) 0 ^ +y f y v + Zi z v ) 
4 lAl 


b t 


1 / “ll “21 


lAl 


I Al/ 


*2* 


1 / “ 1 2 “ 22 ' 
2 K ~ +6 2 — 


lAl 


lAl/ 


(1300 

(131a) 
(131b) 
(131c) 
( 1 3 1 d) 
(1 3 1 e) 
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where 


(1310 


V ~ E p~ \ ^ + y t +z ^ + \ hT A 1 b 


4 lAl " &V\» - x v y i )2 + <Vr, - V« > 2 + <V„ " -Vf )2 > 0 0 

If the results of theorem 1 are used and (£, r\, t) are substituted for (x, y, t ) in equations (72), 


a * j and a* 2 must be of the forms 


n * — 

U \\ 


T^ON^v)- T 2 (t)N { ( 17 ) 


/7 V — 

u 22 


M x (X)N 2 {ri)- M 2 {ON l (t?) 


a *2 =0 


(134a) 


h* = 0 


(134b) 


b* = 0 


It follows from equations (131), (132), and (134) that 


-2^2 = -2«2i = + V* +z f z . = 0 


(134c) 


(135a) 


b x =0 


(135b) 


b 2 = 0 


(135c) 


^ = V 2 (x2+ ^ +z2 > 


From equations (130) and (135), it is seen that 


z i =^+^+Vj =0 


(135d) 


(136a) 


«2 = * f * T , + ¥, + ¥, =0 < 136b ) 

Consequently, to use the results of theorem 1 , the equations of transformation (128) must satisfy 
the conditions 

x t x +y t y + z t z =0 (137a) 

? !) t 1 


x r v +y f y +z t z =0 


(137b) 


x t x + v, v + z t z =0 

t TJ / t-' Tj r TJ 


(137c) 


x 2 +y 2 + z 2 

V y r\ tj 


= a * = 
“ll 


T N - T N 

1 1 JV 2 i 2 iv l 

M^N 2 - m 2 n x 


(137d) 


30 



-T X M 2 + T 2 M X 


(137e) 


x 2 + v 2 + z 2 
y i z ( 

~ a 22 

8 lAl M l N 2 -M 2 N x 

where 

4IAI =(YrVf) 2 +( VrVf ) 2 +(y t z n~ y n z t )2 (138) 

and the potential energy E must be of the form 

E p - \ (x2 t + y 2 + z?)=V= 7*0(0- a* n M 0 &)- a* 2 N 0 ( V ) (139) 

The conditions given by equations (137) to (139) are necessary and sufficient for the Hamilton- 
Jacobi equation describing the motion of a particle on a moving surface to be separable. 

It is obvious that the set of equations (137) and (138) is highly nonlinear and that to use the 
results of theorem 1, one must be able to solve this set of coupled nonlinear first-order partial differen- 
tial equations; actually, at the beginning it is not even known whether any solutions exist. First, let us 
prove the following lemma. 


Lemma 1: Given a Euclidean three-dimensional space E and a set of Cartesian coordinates 
(x, y, z), let (£, t?, t) be new coordinates defined by the transformation 


x = *(£, v, 0 

(128a) 

y - y(£, v, t ) 

(128b) 

z = z(£, TJ, t ) 

(128c) 

Suppose that these equations can be inverted to yield 


t = £0 c,y,z) 

(140a) 

V = r](x,y,z) 

(140b) 

t = t(x,y,z ) 

(140c) 

Then the set of equations 


x t x + y,y + z f z = 0 

f v J v n f r; 

(137a) 

XfX* + y t y % + z f z s = 0 

(137b) 

x,x +y.y + z.z -0 

t T 7 ' t y T\ t T) 

(137c) 

holds if and only if the following set of equations holds 


+f j) +£t? -0 
% x 'x V 'y s z ‘r 

(141a) 

o 

ii 

n 

AW 

+ 

AW^ 

+ 

AW*' 

(141b) 

Wx+tylty +t z%=° 

(141c) 


31 



Proof: Consider the line element 

ds 2 = g.. dx l dx 1 i, j = 1,2,3 

where 


(142) 


x l = y 


x 3 = z 


Note that this may also be written as 


ds 2 = g kl dx k dx 1 k,l= 1,2,3 

where 

(143) 

li 

IK 


X 2 = 1? 


X 3 = t 


Because 


dx‘ dx 1 


Ski ~ Si/ 

dx k dx 1 

(144) 

and 


g.. = 6.. 

°‘i ij 

therefore, 

(145) 

0x dx 3 y d y dz dz 

gl2 ~d^d^ + did^ + d^d^ 

(146a) 

dx dx dy dy dz dz 

gl3 ~ d% dt dj; dt + d£ dt 

(146b) 

dx 3x dy dy dz dz 

gl3 3tj dt 3q dt 3t? dt 

(146c) 

Hence, from equations (137) and (146), it is seen that 


n 

o 

«*«». 

If 

(147) 

so that the coordinate system (£, rj, t ) is orthogonal. That is, the three families of the surfaces given by 

%(x,y,z) = c l 

(148a) 

riix, y, z ) = c 2 

(148b) 

t(x,y, z) = c 3 

(148c) 
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where c { , c 2 , and c 3 are arbitrary constants, are mutually orthogonal. Hence, it follows that 

V£-\7t? = 0 


(149a) 


V^Vt = 0 


(149b) 


Vrj’Vt = 0 


(149c) 


which are precisely equations (141). Conversely, if equations (141) hold, then the three families of 
surfaces given by equations (148) are mutually orthogonal. Consequently, equation (147) holds. Then, 
from equations (144), (145), and (147), it follows that equations (137) hold. This completes the proof 
of the lemma. 

From this lemma, it is seen that to use the results of theorem 1, the equations of transfor- 
mation (128) must be such that one must construct £(x, y, z) = constant and r?(x, y, z ) = constant 
surfaces in the Euclidean three-dimensional space E such that they form a mutually orthogonal set with 
the /(x, y, z) = constant surfaces obtained from the constraint (127). Moreover, these three sets of 
surfaces must also satisfy equations (137d) and (137e). Even though the statement of the conditions 
given by equations (137) has been simplified, the existence of such transformations remains to be 
proved and will be done with the following lemma. 

Lemma 2: Consider a Euclidean three-dimensional space and a set of Stackel coordinates 
(ju, v, gj), which are therefore orthogonal, and let the line element be denoted by 


= 1 g a 


.... dx dx 


(150) 


where 


x - n 


2 

X - V 


3 _ 

X = CO 


Next, consider a particle moving under the action of external forces in the Euclidean three-dimensional 
space but constrained on a moving surface on which the coordinates are (£, r?). If the surface is allowed 
to move in such a way that it assumes the forms of any one set of the Stackel surfaces, and if the 
constant £ and r\ surfaces correspond to the two remaining sets of Stackel surfaces; e.g., 

M = $ (151a) 

v = t\ (151b) 


co =/(0 


(151c) 


where /(/) is a monotonic function (this condition may be relaxed by considering intervals (/j , t 2 ), 
/ 3 ), etc., if /'(/) is not monotonic), then the Stackel coordinate system satisfies the five conditions 
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in equations (137) if and only if the g 33 component of the contravariant metric tensor is not a func- 
tion of co. 


Proof:. For a particle moving in a Euclidean three-dimensional space with a coordinate system 
(n, v, to), it is well known that the coefficients a*, in the Hamilton-Jacobi equation are related to the 
contravariant metric tensor g 27 by the equation 

(152) 


This follows by noting that the kinetic energy E k for a particle of unit mass can also be written as 



1 dx' dx ! 
2 gi > ~dt ~dt 


(153) 


and then making use of equations (21), (30), (36), and (37). (For the case of dynamical problems with 
n degrees of freedom, and not merely those involving the motion of a particle in a three-dimensional 
space, it is also possible to consider an n space and associate with it a metric defined by equation (152).) 
Also, it is noted that for a particle moving in a Euclidean three-dimensional space, with a Stackel 
coordinate system (p, v, co), the coefficients afj in the Hamilton-Jacobi equation may be obtained as a 
special case of the results of theorem 2, n - 3, by letting T 3 = 1 and T { = T 2 = ■ • • = T n = 0 in equa- 
tion (124). Then, 


a* i (M, v,6j)= — 
a\ 2 (M, v, w)= ^ 

i 

a* 3 0 i, v, w)= ~ 

where 



(Gj* 2 

- Ri q 2 ) 

(154a) 

( Ri p 2 

~ p i 

R 2 ) 

(154b) 

PiQi 

- Q\ 


(154c) 

Qi 





*2 


(155) 

Q 3 

R 3 




P t , Q t , and R ; are functions of fx, v, and tu, respectively. Next, consider the motion of a particle con- 
strained on a moving surface in the Euclidean three-dimensional space and choose (£, 77 , t ) as described 
before by the Stackel surfaces and expressed in equations (151). For this, the Hamilton-Jacobi equa- 
tion is given by 


(dS\ 2 

a* n ^,v,m) y—J 


+ a* 2 [%, 


\2 dS 

j + v[S,v,f(. 01 + — =° 


(156) 
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The conditions in equations (137a), (137b), and (137c) are automatically satisfied because all the 
Stackel coordinate systems are orthogonal. Consequently, it remains to prove the following statement: 
The Stackel coordinate system has the property that it results in having the coefficients [£, i?,/(t)] 
and a* 2 It, P,/( 01 in equation (156) satisfy the two remaining conditions (137d) and (137e) if and 
only if g 33 is independent of co. To prove that it is necessary, it is observed that [£, i?,/(t)] and 
a* 2 [£, rj,/(t)] have the form 

T 1 UW 2 (V)~ T 2 (t)N { ( V ) 

a* n U, V, fit)] = (1 57a) 

M^N^ir,) 

-T l it)M 2 a) + T 2 it)M l (^ 

a| 2 [^T7,/(0] = (157b) 

M l i^)N 2 iri)-M 2 iON l i V ) 

From studying equations (151), (154), (155), and (157) and noting that 

A = R x iP 2 Q 2 - Q 2 P 2 )~ R 2 {P x Qs - Q 1 P 3 ) + R 3 (P l Q 2 - Q { P 2 ) (158) 

it is concluded that /?.( co) must be a constant for i = 1 , 2, or 3 and the remaining terms vanish in equa- 
tion (158). Hence, from equations (154) it is seen that a% 3 (ji, v, co) cannot contain co explicitly. That 
is, g 33 is independent of co. Conversely, to prove that it is sufficient, it is observed that if g 33 is 
independent of co, then a* 3 (p, v, co) does not contain co explicitly. Without going into the details of 
Stackel coordinate systems, it cannot be concluded at this point, observing from equations (154), 
(155), and (158), that 7? ( .(oo) is constant for / = 1, 2, or 3 and that the remaining terms vanish in A so 
as to arrive at equations (157). However, if cases are constructed in the way described by equations 
(151) for all the 1 1 Stackel coordinate systems (see section entitled “Results”), it is observed that if 
a| 3 is not a function of co, then equations (157) are satisfied. This completes the proof of the lemma. 

Next, it will be proved that the preceding admissible cases exhaust all possibilities for all coordi- 
nate systems besides those of Stackel. That is, no other coordinate system can give rise to a choice of 
(£, V, 0 such that the five conditions in equations (1 37) are satisfied. This leads us to the following 
lemma. 

Lemma 3: Consider a Euclidean three-dimensional space and a set of orthogonal coordinates 
( fj . , v, oj). Next, consider a particle moving under the action of external forces in the Euclidean three- 
dimensional space but constrained on a moving surface on which the coordinates are (£, t?). Let the 
surface move in such a way that it assumes the forms of a given set of the coordinate surfaces and also 
let the constant | and ry surfaces correspond to the two remaining sets; e.g. 


M= 1 

(151a) 

v = rj 

(151b) 

' \ 

< 

n 

3 

(151c) 


where fit) is a monotonic function. (This condition may be relaxed by considering intervals (fj , t 2 ), 
( t 2 , t 3 ), etc., if fit ) is not monotonic.) If, by this construction, the five conditions in equations 
(137) are satisfied, then the given coordinate system in, v, co) must be of the Stackel type. 


35 



Proof: By starting with an orthogonal system, the three conditions in equations (137a), (137b), 
and (137c) are automatically satisfied. If one now has [£, i?,/(t)] and a | 2 [£, q,/(t)] given by 
equations (157), then it is obvious in view of equations (151 ) that P ( { ju), Qj(v), and R { ( cu) for 
i = 1,2,3 can be chosen such that u^fp, v , co), a| 2 (p, v, cu), and a* 3 (p, v, co) will have the forms 
given by equations (154). From these results and equation ( 1 52), the contravariant metric tensor g' ! 
may be obtained. Now, by the theorem of Eisenhart (ref. 10), it is known that the 1 1 Stackel coordi- 
nate systems are the only real ones in Euclidean three-dimensional space that haveg^ of this form. 

This completes the proof of the lemma. 

The discussion in lemmas 2 and 3 and a study of the section entitled “Results” lead to the 
following theorem. 

Theorem 3: The Hamilton-Jacobi equation 

o (!|) +a 22 +V(i,r\,t)+ — = 0 ( 159 ) 

having a complete integral of the form 

5 ( 1 , 17 , t) = M(£,a) + N(r},a)- T(t,a) (160) 

describes the motion of a particle constrained on a moving surface in the Euclidean three-dimensional 
space if and only if the moving surface assumes the forms of a set of coordinate surfaces (e.g., cj = /(?)) 
in a coordinate system (jc 1 , x 2 , x 3 ) = (p, v, cu) in which the component of the contravariant metric 
tensor g 33 is independent of oj, and p and v are chosen such that p = £ and v = rj. There are only five of 
these surfaces (see fig. 2 ): 

(1 ) The plane in transverse motion 

(2) The plane in rotation about a longitudinal 
axis 

(3) The cylinder in radial motion 

(4) The spherical surface in radial motion 

(5) The cone whose angle is variable 

In these cases, it is not necessary for the motion of 
the surface to be uniform. 

Remarks: A study of the results in the next sec- 
tion reveals that in some cases it is possible to choose 
more than one set of coordinates (£, 77 ) on the surface 
to have the desired form required for separation of 
variables. Every problem involving a particle con- 
strained on a moving surface, satisfying only some 
very general conditions (eq. 1 27), has a Hamilton- 
Jacobi equation that can be expressed in the form of 
equation (1 59) by a suitable choice of an orthogonal 
set of coordinates (£, 77 ) on the moving surface. 



(1) PLANE (2) PLANE 



(3) CYLINDER (4) SPHERE (5) CONE 

Figure 2.— Permissible surfaces. 
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A Time-Dependent Problem 

Consider a particle constrained to move on the surface of a sphere whose radius is r = /(?), where 
f(t) is an arbitrary continuous function of time. It is desirable to investigate in this case the form of 
the coefficients a* x and a| 2 and the potential energy that will permit the separation of the Hamilton- 
Jacobi equation. 

For spherical coordinates ( r , 6, 0), the components of the covariant metric tensor g are given by 

g„= 1 061a) 

g ee =t 2 (161b) 


g 00 =r 2 sin 2 0 (161c) 

The kinetic energy E k of a particle of unit mass is given by equation (153) as 

0 62) 

Because r = f(t), 

E k = + f 2 + / 2 sin 2 6 0 2 ) (163) 


Suppose that b' - 0 in equations (16) and (17). Hence, 
Next, from equations (26) and (27), 


H = 


2fHl) 


P 1 , 


2 f 2 (/) sin 2 d 


pI - \ i Am 2 +E p 


(164) 

(165) 


Consequently, the Hamilton-Jacobi equation is given by 

1 /dS\ 2 1 /dS\ 2 i 95 

— — + — )--f 2 +E n + — = 0 (166) 

2 f 2 \d0 / 2 f 2 sin 2 d \90/ ^ 3 1 

which may also be obtained more directly from equations (1 52) and (156). Note that the coefficients 
of (9 5/90) 2 and (35/90) 2 are of the forms required in equations (137), which are 

T x {t)N 2 ^)- T 2 (/W 1 (0) 

a* 0 = i (167a) 

M 1 (0)7V 2 (0)-M 2 (0)/V 1 (0) 


This can be seen by choosing 


-7’ 1 (/)M 2 (0) + r 2 (/)A/ 1 (0) 
M 1 (0F 2 (0)-M 2 O9W 1 (0) 
1 

T \ = 

2/ 2 (0 


( 1 67;b) 

(168a) 
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Mj = 1 (168b) 

TVj = 0 (168c) 

r 2 =0 (169a) 

1 

M,= (169b) 

sin 2 6 

N 2 = 1 (169c) 

Furthermore, E p is chosen to satisfy equation (139) so that 

e „ - \iAm 2 - r 0 w- «;,*,(»)-«») (no) 


where TAt), MA6), and N n (<p) are arbitrary functions of their respective arguments. Consequently, a 


complete integral of the form 

S(d,4>, r,a l ,a 2 ) = M(d,a l ,a 2 ] 

is obtained for the Hamilton-Jacobi equation (166). 
are determined from the equations 



+ N(<t>, aj , Q! 2 ) - T{t, , a 2 ) (171) 

From equation (11), the functions M, N, and T 

a 2 

a. (172a) 

sin 2 9 

a 2 (172b) 


RESULTS 


dT a 2 

= T 0 + (172c) 

dt 2f 2 (t) 


a 


* 

ll 


Based on the discussion in the previous section, results are presented summarizing the forms of 
a* 2 , and E p , where 

E p = V - (x 2 + y] + z 2 ) 


for which one has a separable Hamilton-Jacobi equation. 


In the following subsections describing each of the 1 1 Stackel coordinate systems, h f for i= 1,2,3 
denotes the scale factors, which are related to the components of the covariant metric tensor g 
defined by 


ds 2 = 



(173) 
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and A k for k = 1, 2, . . . , 1 1 denotes the 1 1 Stackel determinants defined by equations (154) and 
(155). These determinants are the transpose of those found in Morse and Feshbach (ref. 1 7), from 
which have been taken geometrical representations of the various coordinate systems. 

Rectangular Coordinates (1) 

h x = h 2 = h 3 = 1 

y = ^2 
Z = h 

Unconstrained Particle 


z 



*1 = 


1 0 0 

-1 1 0 

-1 0 1 


= 1 


= /7* — /7* = 1 

a n a 22 a 33 1 


V=X 0 (x)+Y 0 (y)+Z 0 (z) 

X Q (x) Y Q (y) Z Q (z) 

= + + 


fit 


h\ 


h\ 


Constrained Particle on a Moving Surface 


(1) If 


X=fy (f) 

y = y 


then 


(2) If 


z = z 


/7* = n * =: 1 

a 22 u 33 1 


E p ~ f\ + Tq( 0 + Y Q (y) + Z Q (z) 


x = x 
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y = f 2 (0 


then 


(3) If 


z - z 


a* -a* =1 
“ll u 33 1 


E =f 2 + T 0 (t) + X 0 (x) + Z Q (z) 


X = X 


y = y 

z-/ 3 (0 

then 

n * = n * — 1 

a ll U 22 1 

E p = f]U) + T Q (t) + X 0 (x) + Y 0 (y) 

Cylindrical Coordinates (2) 

h 2 = h\ = 1 

* 2 2 = «? 

X = COS ^ 

^ sin £ 2 
z = £ 3 

*1=' 

^3 =Z 

?! -(x 2 +y 2 ) 112 



tan £ 2 


7 

z 


£3 = Z 
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Unconstrained Particle 


^2 = 


1 0 0 


1 0 


-1 0 1 


= 1 


a 


* 

li 


= 1 


1 



V-X Q {r)+ — “ — +Z Q (z ) 


X 0 (r) Y o m 

h\ h\ 


Z 0 (z) 


Constrained Particle on Moving Surface 


(1) If 


then 


(2) If 




e = e 


z = z 


*22 




Ep f\ + T Q {t) + a^ 2 Y q (.®) + a*'iZ Q { z ') 


r - r 


e =/ 2 (/> 

z = z 
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then 


(3) If 


*11 = * 3*3 = 1 


E p =r 2 f 2 2 +T 0 {t) + X 0 (r) + Z 0 (z) 


Spherical Coordinates (3) 


(I * ~ 1 

a n 1 


n ^ — 

a 22 


Z=f 3 0) 


K =/] +r 0 (O + < 1 X 0 W+^ 2 l' 0 W) 



*5 = i 


*2 = «| 

^3 = I? sin 2 

x = sin £ 2 cos £ 3 

^ sin £ 2 sin £ 3 

z = COS £ 2 


*1=' 

£, =<? 


^ = r = (x 2 +y 2 + z 2 ) 1/2 


cos = — 
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Unconstrained Particle 


tan £ 3 


y 

x 


*3 = 


-1 


sin 2 6 


= 1 


7 * — 

^11 


7 H 8 = 

*22 


7 ^ — 

*33 


r 2 sin 2 0 


F=X 0 (r) + -r o (0) + 


/- 2 sin 2 0 


Z n ( 0 ) 


x 0 (r) r o (0) z o (0) 

+ + 


h \ 


Constrained Particle on Moving Surface 

(1) If 


h 2 

n 2 


hi 


r ~f \ (0 
0 = 0 


then 


4> = <p 


a* = — 
22 c2 


1 

A 
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1 


(2) If 


then 



f\ sin 2 0 


^=/?^ 0 W +<I .r oW) + a ,, w 


r~ r 


0s f 2 (O 


0 = 0 



1 


(3) If 





r = r 


then 


0 = 0 


7* =: 1 

l ll 1 



*=f 3 (0 


E P ~ r 2 sin 2 0 /'2 

Elliptic Cylindrical Coordinates (4) 


+ 7 ’oW + af I X o (,) + fl|2ro(0) 



cos 2| 2 ) 
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where d is a constant 


x - d cosh cos £ 2 


y = d sinh ^ sin £ 2 
Z = *3 


where 


Unconstrained Particle 


$ 2 = ^ 

^3 = Z 

sinh 2 = /I +5 
sin 2 £ 2 = -/4 + 5 

^3 ~ z 
1 

i4 = (x 2 +y 2 - d 2 ) 

2d 2 


B = 



+ 


y 2 V 12 

d 2 / 




d 2 

— cosh 2£ 


d 2 

— cos 2 t? 
2 


- 1 


d 2 

= — (cosh 2£ - cos 27?) 


2 

n * = n * = 

11 22 Pr , 

d~ (cosh 2£ - cos 2r?) 


7 * — 

*33 



v= 


2X 0 (Z) 


2 Y q {t}) 


d 2 { cosh 2£ - cos 2r\) d 2 (cosh 2£ - cos 2 tj) 0 ^ 


x 0 (H) r 0 (n ) z 0 (z) 

h 2 + a| + h 2 

Constrained Particle on Moving Surface 

(1) If 

*1 =/,W 

£ 2 

^3 =Z 

then not separable. 

(2) If 

^3 = Z 

then not separable. 

(3) If 

=/ 3 w 

then 

2 

a l 1 = fl 22 = 

d 2 (cosh 2£ - cos 277) 

£ ,=/3 + ro+«f,[-r„tt)+r 0 w] 

Parabolic Cylindrical Coordinates (5) 
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> 


Unconstrained Particle 


*- {«?-{!) 
y = z 1^2 

Z = S 3 

= x + (y 2 - x 2 ) 1/2 
£ 2 = -x + O 2 - x 2 ) 1/2 
^3 ~ Z 


*5 = 


*1 

-1 



fl ll = fl 22 


1 



! 33 = 1 

w 

W 


F= 

+ 

+Z Q (z) 



W 

W 

3>fri ^ 


+ + 


*} h l h \ 


*o«i> 


2(> 2 - x 2 ) 1/2 


20 2 - x 2 ) 1/2 


+ Z 0 (z) 


47 



Constrained Particle on Moving Surface 


(1) If 

%2 = %2 
£3 = Z 

then not separable. 

(2) If 

£ 2 =/ 2 (0 

£3 = Z 

then not separable. 

(3) If 

= ^2 


then 


Conical Coordinates (6) 


? 3 = AW 

1 

Z *1 “ a 22 = T 7 

*I + *2 

^=/?(0 + 7’ 0 +«f 1 iJr 0 (£i) +ir o(£ 2 )] 



h\=h\ 

= £ 2 [A : 2 cn 2 (|j , A) + A ' 2 cn 2 (| 2 , A')] 
A 2 = 1 

where A and A' are constants of the elliptic 
functions and 

A 2 + A ' 2 = 1 
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* = £ 3 dn (£j , A) sn (£ 2 , A') 
.A = £3 sn (£j , A) dn (£ 2 , A') 


Because 


z = £ 3 cn (^,k) cn (% 2 ,k') 


where 


Unconstrained Particle 


sn 2 - cn 2 = 1 
dn 2 + k! sn 2 = 1 


dn 2 = k 2 cn 2 = k 
k 2 cn 2 (£p k) = A + B 
k ' 2 cn 2 (£ 2 > k') = -A + B 


1 

A = — [A 2 * 2 - At' 2 3> 2 + ( k 2 - k ' 2 )z 2 ] 
2r 2 


B = 



k 2 k' 2 z 2 \ 1/2 



A 6 = 


0 0 

k 2 cn 2 (| ls A) k ' 2 cn 2 (£ 2 , A') 
-1 1 


t 2 

-1 

01 


= £3 [A 2 cn 2 , k) + k ' 2 cn 2 (£ 2 , k')] 

1 


//* =/ 7 * = 

W 11 a 22 


£3 [A 2 cn 2 (£j , A) + k ' 2 cn 2 (£ 2 , A:')] 


= 1 

“33 1 




^0 ^0 

- — + — + z o 

h\ h 2 



Constrained Particle on Moving Surface 


(1) If 

€i=/i 
%2 ~ %2 
£3 = ^3 

then not separable. 

(2) If 

^2 =/ 2 W 

^3 = ^3 

then not separable. 

(3) If 

*i=«i 
£2 = ^2 

€3 “AW 

then 

1 

^ * = Q* 

/ 2 [Ac 2 cn 2 (£, , A:) + k ' 2 cn 2 (| 2 , &')] 

E p =/](0 + ^(0 + a*, [* 0 ($, ) + r 0 (| 2 )] 

Parabolic Coordinates (7) 
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*V 


where 


e 2 =r- z 

y 

tan £, = — 

J v 


r = (x 2 + y 2 + z 2 ) 1,)2 


Unconstrained Particle 


$ 


2 

1 


& 


0 


1 1 



-1 1 0 



2 

2 


* = a*- — 


‘11 


22 


* l + f 


a 


* 

33 


1 

$8 


v = 


*0«1> r 0«2> Z 0^3> 

+ + 


*1 + *2 




fc 2 >2 
?1 ? 2 


^ + ^o + fo 

h 2 hj h l 


Constrained Particle on Moving Surface 


(1) If 


«1 = /.w 


l 2 ^2 

? 3 = £3 
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then not separable. 

(2) If 

i 2 = / 2 w 

^3 = ^3 

then not separable. 

(3) If 

%2 = 

h “/ 3 « 

then 


— /j'i> “ 

*11 a 22 


fl + fl 


= ?h^(0 + r 0 +a«, [*<,({,)+ r 0 « 2 )] 

Prolate Spheroidal Coordinates (8) 



h\= h\= d 2 (sinh 2 % } + sin 2 £ 2 ) 
h\ = d 2 sinh 2 ^ sin 2 £ 2 
x-d sinh sin £ 2 cos £ 3 

^ sinh sin £ 2 sin £ 3 

z = d cosh cos £ 2 


sinh 2 £j = zl + B 


sin 2 % 2 - B - A 


where 


tan £ 3 


y_ 

x 


1 

A = (r 2 - d 2 ) 

2d 2 
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B = 




1/2 


r 2 = x 2 +y 2 + z 2 


Unconstrained Particle 


d 2 sinh 2 

d 2 sin 2 ^ 2 

0 

1 

1 

-1 

sinh 2 

sin 2 

-1 

1 

0 


- d 2 (sinh 2 + sin 2 £ 2 ) 


1 

= (l * ‘ = — — 

<i 2 (sinh 2 + sin 2 £ 2 ) 

1 

d 2 sinh 2 sin 2 !~2 

r 0«2> Z 0« 3 ) 

— + + 

d 2 (sinh 2 £j + sin 2 £ 2 ) c/ 2 (sinh 2 ^ + sin 2 £, ) t/ 2 sinh 2 sin 2 £ 2 

_ z o + ^o + fo 

*1 + h l + h \ 


Constrained Particle on Moving Surface 


(1) If 



then not separable. 


tj'fj t 3 =t: 

(2) If 



then not separable. 


^2 ~ f 2^} ^3 = £; 
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(3) If 






£3 =/ 3 <') 


“ll “22 


d 2 (sinh 2 + sin 2 £ 2 ) 

£ = d 2 sinh 2 £j sin 2 £ 2 / 2 + r o (t) + a*j [3f 0 (£j ) + F Q (£ 2 )] 


Oblate Spheroidal Coordinates (9) 

z 



h 2 = h 2 = d 2 {c osh 2 £j - sin 2 £ 2 ) 
ft 2 = d 2 cosh 2 sin 2 £ 2 
x = d cosh sin £ 2 cos £ 3 
y = rf cosh £, sin £ 2 sin £ 3 
z-d sinh £j cos £ 2 

cosh 2 = /4 + B 
cos 2 i 2 = - A + B 

. , _ y 


where 


Unconstrained Particle 


A= — O 2 - d 2 ) 
2d 2 


2 \ 1/2 


B= \A Z + 


d 2 cosh 2 

d 2 sin 2 % 2 

1 

1 

cosh 2 

sin 2 jf 2 

-1 

1 


1 = d 2 (cosh 2 + sin 2 £ 2 ) 
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1 


fljl ° 22 d 2 (cosh 2 + sin 2 t 2 ) 
1 


d 2 cosh 2 ^ sin 2 £ 2 

* 0^1 > + ^2 > Z 0 ^ 


d 2 (cosh 2 + sin 2 £ 2 ) d 2 cosh 2 £j 

Constrained Particle on Moving Surface 

0) If 

^2 = ^2 
^3 " ^3 

then not separable. 

(2) If 

^2 = A 

^3 = ^3 

then not separable. 

(3) If 

^2 = ^2 

S 3 “/ 3 W 

then 

1 

Q * = (l * = 

d 2 (cosh 2 + sin 2 £ 2 ) 

E p = cl 2 cosh 2 £j sin 2 £ 2 / 2 + 7’ 0 (f) + af 1 [A^Gf. 


) 

sin 2 £ 2 


l ) + Z o^2 ^ 
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Ellipsoidal Coordinates (10) 



«; - nix # - 1 \) 
(g - gxg ~ g) 

-a 2 )(l| - £ 2 ) 

gi - ilxtl - s 2 2 ) 

(£ 2 " a 2 )(£ 2 - b 2 ) 


^ = 


"(I 2 - 6 2 )(|2 - b 2 )(%\ - b 2 ) 
b 2 (b 2 -a 2 ) 


1/2 


ab 


Unconstrained Particle 


^10 


-'a 3 


g-a 2 


&~a 2 


(£ 2 - b 2 )(a 2 - b 2 ) (X 2 - b 2 ){a 2 - b 2 ) (| 2 - b 2 ){a 2 - b 2 ) 

(*j ~ €|)«j ~ Kg - S?) 

(£ 2 ~ a 2 )(€l " a 2 )^ 2 - fl 2 )(| 2 - b 2 )(X\ - b 2 ){% 2 - b 2 ) 


(Z\-a 2 )(Z 2 -b 2 ) 

(t| - €f )«i - €j) 
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-Gl-cPxtl-b 2 ) 

tt?-« 2 X*f- ft 2 ) 
ttl- $?X$f- *1) 

F - <2j*j ^ 0 (^i ) + <222 ^ 0 ^ 2 ) + a *3 ^0 ^3 ) 

Constrained Particle on Moving Surface 

(1) If 

Si-A(') 

£2 = ^2 
£ 3 = ^3 

then not separable. 

(2) If 

*i=*i 

S 2 =/ 2 (0 

£3 = ^3 

then not separable. 

(3) If 

*i“*i 
£2 = 

*3 “W 

then not separable. 

Paraboloidal Coordinates (11) 

, , - €iX€? - €f) 

h \ = £j 

(%\-a 2 X%\-b 2 ) 

2 _ gljjixgjS) 
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h\ =& 


x = 


y = 


t «S - gjxtj - tj) 

! (^3 - - ^ 2 ) 

(*? - fl 2 )«l ' A 2 )« 2 3 - A 2 )' 
a 2 - b 2 

\k\ - b 2 M\ - b 2 )tt 2 - b 2 ) 


1/2 


1/2 


b 2 - a 2 


Unconstrained Particle 


A ll = 


* = 2 «l + « 2 + «3 


Z> 2 ) 


«? 

*? 


^2 


s 2 - * 2 
*1 


£ -« 2 


£ 


fc2 

*3 


3 


t 2 - / j 2 


fc2 

*3 


(^ 2 - b 2 )(a 2 - b 2 ) 




b 2 ){a 2 - b 2 ) (g-b 2 )(a 2 -b 2 ) 


(S 2 3 - €lX«J - €?X «1 - 5 ? ) 


= *1*2*3 


‘11 


22 


(*J - « 2 )(| 2 2 - « 2 )« 2 3 - a 2 )«? - 6 2 )(S 2 - 6 2 )(t 2 - /> 2 ) 
1 (i\-a 2 )& 2 -b 2 ) 

*? ai - * lx *? - « l ) 

1 «|-fl 2 )«|-6 2 ) 

1 ~a 2 ){i 2 -b 2 ) 


a h = 


H « 1 -€? X « 2 3 -€|) 

P = af ! , €2 ’ *3 )* 0 tti ) + *2*2 (*! ^2 . $3 )7 2 (*2 ) + a 3*3 ^1 ’ *2 > *3 > Z 0 ($3 
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Constrained Particle on Moving Surface 


(1) If 

=/ 1 W 
^2 = ^2 
^3 = ^3 

then not separable. 

(2) If 

«2 a AW 

* 3 = *3 

then not separable. 

(3) If 

^2 = ^2 

^3 = / 3 « 

then not separable. 

DISCUSSION 

Application of Present Results to Other Classes of Dynamical Problems 

In a previous section, the results of theorem 1 for the time-dependent two-dimensional Hamilton- 
Jacobi equation were applied to the solution of the class of problems describing the motion of a 
particle constrained on a moving surface in a Euclidean three-dimensional space. 

Theorem 2, n - 3, may similarly be applied to study problems comprising two moving particles 
subject to three constraints. For example, consider the following class of problems: Let/ > 1 andF 2 
respectively denote two particles with coordinates (x { ,y l> z l ) and (x 2 ,y 2 ,z 2 ). Suppose these 
particles are constrained on the same moving surface and also the distance between them is constant. 
Then, 


f(x l ,y l ,z l , 0 = 0 
f(x 2 ,y 2 ,z 2 , t) = 0 
(*1 ~ x 2 ) 2 +(y 1 - y 2 ) 2 + (z, - z 2 ) 2 = d 2 


(174a) 

(174b) 

(174c) 
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and the constraint may be used to reduce the given problem to one involving only three independent 
generalized coordinates (£, r?, £)• Using the conditions given in theorem 2, n = 3, a general class of 
moving surfaces may be found that leads to problems solvable by the separation of variables in the 
Hamilton-Jacobi equation. 

. In this way, the comparatively wide classes of problems that can be solved by the results of 
theorem 2 (which deals with n degrees of freedom) may be studied. There is no confinement to prob 
lems in Newtonian mechanics because the results are purely mathematical. Thus, problems arising in 
relativity and even any other non-Newtonian form of mechanics may be considered as long as one has 
a Hamilton-Jacobi type of equation. 


More General Forms of the Hamilton-Jacobi Equation 


The Hamilton-Jacobi equation studied in this work is diagonal and also does not have first- 
degree terms. Thus, one may consider the general diagonal Hamilton-Jacobi equation 


(60) 


/ 35 \ 2 

^ bS 

bS 

I — I + 

> b * — 

+ V+ — = 0 

W/ 

fk a?, 

bt 


;= l 

or the most general Hamilton-Jacobi equation 

V V * 35 

2 L 2 -, v 3 q. bq. Z-> ' 3 q. 

(=1 /= 1 » = 1 

and seek complete integrals of the form 


3 S 

+ > b? — + V + — = 0 


S b < 


3 S 

Jt 


n 

SXq, f;a) = ^ X j (q i ; a) - T{t\ a) 

i= l 


(54) 


(ID 


where X i and T are arbitrary functions of q { and t. These have already been done by Chan. 1 


Other Forms of Separation 

Ever since the time of Jacobi, complete integrals of the Hamilton-Jacobi equation have always 
been sought in the form of equation (11). In fact, the term “separation of variables” is defined spe- 
cifically with reference to this form. It has long been known that many dynamical problems have 
Hamilton-Jacobi equations that do not possess complete integrals of this form. However, essentially 
all the problems illustrated in texts that are solvable by separating variables of the Hamilton-Jacobi 
equation can also be solved without the sophistication developed in the theory of Hamilton and 
Jacobi. For example, one may solve the problem of the linear harmonic oscillator in many steps 
fewer and more simply by starting from Newton’s equations than by starting from the Hamilton- 
Jacobi equation. But this lack of usefulness of the Hamilton-Jacobi theory is only superficial. 


1 F. K. Chan: 


“Classical and Non-Classical Separation of the Hamilton-Jacobi Equation,” in preparation. 
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In accordance with terminology introduced by Chan, 2 separation of variables yielding complete 
integrals of the form <pf equation (11) shall be referred to as classical separation and separation of 
variables yielding complete integrals not of the form of equation (1 1) shall be referred to as nonclassi- 
cal separation. Thus, as an example, he has considered the general Hamilton-Jacobi equation (54) and 
sought a complete integral of the form 

n 

S(q, t',a)~ ^ T(^;a) - T(t\ a) + F(q, /) (175) 

i= l 

where F(q, t ) is an arbitrary function specified a priori. 

Disadvantages of Specifying Forms of Complete Integrals 

Even though the form of equation (175) is more general than the form of equation (1 1 ), it turns 
out that it still does riot yield the solutions of many interesting dynamical problems. In fact, this is 
the main disadvantage in the whole outlook of seeking complete integrals of some specified form. For, 
even after one has obtained necessary and sufficient conditions for such separation, one has then to 
find physical problems that satisfy these conditions. Whereas, in life, one is usually given a definite 
physical problem to which one has to find the answer. 

Thus, after considering many forms of complete integrals, Chan 3 is still unable to handle problems 
involving forces resulting from rotating sources. Specifically, there is the time-honored circular 
restricted planar three-body problem described by the Hamilton-Jacobi equation 

1 

+ - 
2 

where 

P? = ($ - a) 2 + V 2 (177a) 

f>\ =(£ + 6) 2 +V 2 (177b) 

for which a complete integral involving a finite number of terms is sought. However, he has been able 
to obtain solutions involving a finite number of terms for problems with rotating sources if the poten- 
tial is given by 

n i 

v = Y -kfKS-a;) 2 +07- bfl (178) 

i=i 

which is of the type described by Hooke’s law. a 

Separation of Variables in Higher Order Partial Differential Equations 

The technique of separation of variables has long been applied to the solution of the usual linear 
second-order partial differential equations of mathematical physics: 

2 F, K. Chan: “Classical and Non-Classical Separation of the Hamilton-Jacobi Equation,” in preparation. 

3 F. K.Chan: “Reformulation of Problems With Rotating Force Fields.” Submitted to Celestial Mechanics. 



/a w\ 2 

a w 

a w /gm 1 

co£ 

gm 2 \ 


+ cor? 

+ — H 

(176) 

A 37?/ 

9£ 

\P! 

p 2 / 
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Heat equation: 

, 1 d u 

V 2 u = 

a 2 dt 

(179) 

Wave equation: 

, 1 d 2 u 

V 2 u 

c 2 dt 2 

(180) 

Laplace equation: 

V 2 u = 0 

(181) 

Helmholtz equation: 

V 2 u + \u - 0 

(182) 

Klein-Gordon equation: 

Vu-Pw ' S2 “ 

c 2 dt 2 

(183) 

Schrodinger equation: 

- h 2 , T , du 

V 2 u + Vu = ifi — — 

2 M dt 

(184) 

Thus, one seeks a solution of the form 



u(x, t ) = T(t-a) JJ X fa-fit) 

i=\ 

(185) 

or 

u(x) = jj X^x-ia) 
i= 1 

(186) 


and substitutes into the differential equation in question to verify that a solution of this form exists. 

In physical situations, the coordinate system is usually one of the 1 1 Stackel forms and, consequently, 
the question of separability has been answered completely by Robertson (ref. 9) and Eisenhart (ref. 
10 ). 

However, the question of separability for the most general linear second-order partial differential 
equation has not been answered yet in the sense that no necessary and sufficient conditions equivalent 
to Robertson’s for equations involving the Laplacian operator have been published. Thus for the sake 
of convenience, if one considers the forms 


n 


d 2 u 

n 


du 


du 



E 

a.(x, 0 

— 

-E 

h (x, t ) 

— 

+ c(x, t)u - 

— 


(187) 

1=1 


dx { 2 

/=1 


dx i 


dt 



n 


d 2 u 

n 


du 


d 2 u 

du 


E 

i=i 

a (x, t ) 

dx 2 

+ E 

i=i 

h, (x, t ) 

dx i 

+ c(x, t)u -- 

= + d(x, t) 

dt 2 

dt 

(188) 
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" 9 2 m « 

Y! a ,( x ) — - + 


E 6 « (x) 


9m 


+ c(x)m - 0 


bx i 


(189) 


r-f ' 9* 2 . 

1=1 i i=i 

many problems of physical interest are included, such as the three-dimensional heat equation with 
space- and time-dependent specific heat c(x, t) and density p(x, t ) and also the equation for non- 
isotropic, space- and time-dependent thermal conductivity fc.(x, t): 

3 1 9 / 9 u \ 9 u 

(190) 


r- > 3 / 3 U \ 


9m 

"9? 


the one-dimensional wave equation for a vibrating string with space- and time-dependent tension 
p{x, t ), elastic restoring force q(x, t), and density p(x, t ): 

9 2 m 

- p(x, f)« = p(x, t) 


9x 


9m 

p(x, 0 — 

9jc 


(191) 


9t 2 


and the two-dimensional wave equation for a vibrating membrane with space- and time-dependent 
density p(x, t) and nonisotropic, space- and time-dependent tension p ( (x, t ): 

2 ''l - 9 m"| 9 2 m 

P,-(x, 0 — 


9 


9x 


= p(x, t ) 

9 1 2 


(192) 


The problem of separation of variables in the forms of equation (185) or ( 1 86) for the most 
general linear second-order partial differential equation expressed for convenience in equations (187) 
to (189) has also been studied by Chan. 4 


CONCLUSIONS 

Theorem 2 gives necessary and sufficient conditions for the separation of the special diagonal 
Hamilton-Jacobi equation. Theorem 3 gives the complete set of solutions of the dynamical problem 
of a particle constrained on a moving surface solvable by the separation of the special diagonal 
Hamilton-Jacobi equation. There are only five classes of such constraining surfaces: 

(1) The plane in transverse motion 

(2) The plane in rotation about a longitudinal axis 

(3) The cylinder in radial motion 

(4) The spherical surface in radial motion 

(5) The cone whose azimuthal angle is variable 

In these cases, it is not necessary for the motion of the surface to be uniform. 

Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland, July 11, 1972 
311-80-43-01-51 

4 F. K. Chan: “Separation of a General Second-Order Linear Partial Differential Equation of Mathematical Physics,” in preparation. 
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